Electronic structure of the quasi-one-dimensional organic conductor TTF-TCNQ 
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We study the electronic structure of the quasi-one-dimensional organic conductor TTF-TCNQ 
by means of density-functional band theory, Hubbard model calculations, and angle- resolved 
photoelectron spectroscopy (ARPES). The experimental spectra reveal significant quantitative and 
qualitative discrepancies to band theory. We demonstrate that the dispersive behavior as well as the 
temperature-dependence of the spectra can be consistently explained by the finite-energy physics of 
the one-dimensional Hubbard model at metallic doping. The model description can even be made 
quantitative, if one accounts for an enhanced hopping integral at the surface, most likely caused by 
a relaxation of the topmost molecular layer. Within this interpretation the ARPES data provide 
spectroscopic evidence for the existence of spin-charge separation on an energy scale of the conduc- 
tion band width. The failure of the one-dimensional Hubbard model for the low-energy spectral 
behavior is attributed to interchain coupling and the additional effect of electron-phonon interaction. 



PACS numbers: 71.20.Rv, 79.60.Fr, 71.10.Pm 

I. INTRODUCTION 

The electronic structure of one-dimensional (ID) con- 
ductors provides a valuable testing ground for the study 
of the quantum-mechanical many-body problem. On 
the theoretical side there exist various models for ID 
interacting electron systems, which predict highly un- 
usual low-energy excitations due to dynamical decou- 
pling of charge and spin degrees of freedom. As a con- 
sequence, the low-energy paradigmatic Fermi liquid pic- 
ture fails for ID metals and a new generic many-body 
quantum state emerges which is commonly referred to 
as Tomonaga-Luttinger liquid (TLL).^ Experimentally, 
quasi- ID metals are indeed found to display marked de- 
viations from conventional metallic behavior, such as the 
absence of Fermi-Dirac edges in the single-particle excita- 
tion spectra probed by angle-resolved photoelectron spec- 
troscopy (ARPES) iSi^ However, an unambiguous spec- 
troscopic identification of spin-charge separation and the 
existence of low-energy TLL behavior in ID metals is still 
lacking so far. Additional interest in ID electron systems 
arises from the suggestion that their physics may also 



be relevant to the electronic structure of the cuprate- 
based high-temperature superconductors,'* related to the 
recent discovery of charge ordering in these materials 
into narrow metallic ID stripes separated by insulating 
regions. ^'^ Against this background the search for positive 
spectroscopic signatures of unusual electronic correlation 
effects in ID metals remains to be of topical importance. 

In the search for promising realizations of a proto- 
typical (quasi-) ID conductor the organic charge transfer 
salts appear as interesting candidates. Due to the for- 
mation of linear molecule stacks in the crystal structure 
and an electronic charge transfer from cationic to an- 
ionic complexes they display strongly anisotropic metal- 
lic conductivities iSi Photoemission experiments on such 
materials often find unusual spectral behavior like the 
absence of a metallic Fermi edge.'^' However, the lack 
of information on surface quality, especially with re- 
gard to the rapid photon-induced decomposition of or- 
ganic compounds in the vacuum ultraviolet,^ casts se- 
rious doubts to what extent these observations reflect 
intrinsic electronic properties or rather a strongly dis- 
turbed surface. This is further corroborated by the fail- 
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ure of ARPES to detect spectral energy-vs.-inomentum 
dispersions in most charge transfer salts A 
notable exception is TTF-TCNQ (tetrathiafulvalene- 
tetracyanoquinodimethane) , being the first (and so far 
only) organic conductor for which dispersing ID bands 
have been observed by ARPESi^^^ This indicates a well- 
ordered periodic surface structure and thus lends much 
enhanced significance to the observation of a deep pseu- 
dogap around the Fermi energy, which even increases up 
to room temperature. This spectral behavior has recently 
been interpreted as possible indication of a highly un- 
usual normal state in this ID conductor liii^ 

In this paper we present a comprehensive experimental 
and theoretical study of the electronic structure of TTF- 
TCNQ, elaborating on our earlier ARPES results pub- 
lished in R.efilfii. The comparison between experiment 
and band theory reveals significant discrepancies, con- 
cerning both the width of the conduction bands as well 
as their qualitative dispersion. While the band width 
renormalization can be attributed to a molecular surface 
relaxation, the remaining discrepancies indicate a fail- 
ure of the bare band picture. Rather, we are able to 
demonstrate that the TCNQ-derived part of the ARPES 
finite-energy dispersions can be consistently mapped onto 
the electron removal spectrum of the ID Hubbard model 
at finite doping. ^^'^^ The importance of electronic cor- 
relations is further corroborated by a peculiar temper- 
ature dependence of the photoemission spectra. Based 
on these findings the spectral behavior of TTF-TCNQ 
is interpreted as spectroscopic evidence for spin-charge 
separation on an energy scale as large as the conduction 
band width. 



II. PROPERTIES OF TTF-TCNQ 

The monoclinic crystal structure of TTF-TCNQ is 
shown in Fig. ^ The lattice parameters at room tem- 
perature are a = 12.298 A, 6 = 3.819 A, and c = 18.468 
A, the monoclinic angle is (3 = 104.46°>i^ The impor- 
tant structural features are parallel linear stacks of planar 
TTF and TCNQ molecules, respectively, oriented along 
the crystallographic b direction. The 7r-type molecular 
orbitals, extending over the entire size of each molecule, 
overlap with those of the neighboring molecules stacked 
above and below. Maximum covalent bonding is achieved 
by tilting the molecular planes slightly about the a-axis, 
by i9f = 24.5° and i9q = 34.0° for the TTF and TCNQ 
stacks, respectively (cf. Fig. ^J^^ The sign of the tilt 
angle alternates between neighboring stacks, leading to 
the herringbone structure of Fig. 

As covalent bonding occurs only along the stack di- 
rection, the corresponding electronic TTF and TCNQ 
bands are expected to be strongly anisotropic. Charge 
transfer of ^ 0.59 electrons per molecule from TTF to 
TCNQ drives both types of chains metallici^iiS The con- 
ductivity along b is two to three orders of magnitude 
larger than perpendicular to it, making TTF-TCNQ a 



truly quasi-lD metal. Below Tp = 54 K a charge den- 
sity wave (CDW) develops along the b-direction, with 
wave- vector Qcdw = 0.295 b* (= 0.485 A^i). The oc- 
curence of the CDW is accompanied by a metal-insulator 
transition. From the activated behavior of the conduc- 
tivity a Peierls gap of 40 meV has been inferredi^ 
Within mean-field weak coupling theory this translates 
into a transition temperature of Tmf ^ 125 K. Due to 
the dominant effect of fluctuations in ID systems the ac- 
tual Peierls transition is suppressed to about half of this 
value. An additional transverse ordering transition oc- 
curs at 38 Ki^ The observation of diffuse x-ray scattering 
at Q = Akp up to 220 K indicates the importance of 
electronic correlations 4? 



III. BAND STRUCTURE CALCULATION 

The theoretical band structure was studied within the 
standard density functional theory (DFT) approach us- 
ing the generalized gradient approximation (GGA)iiS 
We used the self-consistent full-potential linearized aug- 
mented plane wave (LAPW) method as implemented in 
the WIEN97 code.^o A basis set of about 12500 LAPWs 
and additional "local orbitals" for the 2s (3s) states of C 
and N (S) were employed. This corresponds to a lower 
basis set convergence than desirable, but was limited by 
the available computational resources. Self-consistency 
was achieved using 18 k-points in the irreducible wedge 
of the BZ and a temperature broadening scheme with 5 
mRy. The results are largely consistent with previous 
band calculationSfSLSSiS^ but contain more detailed and 
reliable information due to the more advanced method. 
Very good agreement is found with the recent LDA/GGA 
pseudopotential calculation of Refl23. Our calculations 
have been performed for the experimental room temper- 
ature structure. In addition, we also studied the effect 
of structural distortions as model for a possible surface 
relaxation (see section IVl|) . 

According to the DFT calculation TTF-TCNQ is char- 
acterized by strong intramolecular covalent bonding, 
whereas the interaction between the molecules is predom- 
inantly ionic. Thus, the molecular orbitals are strongly 
localized except along the stacking direction, where small 
but notable covalent intermolecular bonding occurs. As 
seen in Fig. [3 this leads to the formation of two sets 
of quasi-lD conduction band doublets with pronounced 
dispersion along b*, i.e. the FZ line of the Brillouin zone 
(cf. Fig.^. The first one, just below Ep at the F-point 
and unoccupied at Z, is derived from 7r-bonded 2p or- 
bitals of mostly C(6) atoms (in the notation of Ref. fl^ 
and can thus be attributed to the TCNQ stacks. The 
130 meV splitting at F results from a weak interaction 
between the two TCNQ stacks in the unit cell. The other 
conduction band doublet, showing the opposite disper- 
sion and remaining nearly degenerate, is mainly derived 
from the 3p7r orbitals of S(l) and S(2) atoms and there- 
fore associated with the TTF stacks. A small hybridiza- 
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tion gap opens between the respective upper TCNQ and 
TTF bands at the Fermi level. The definite assignment of 
each band to either the TCNQ or the TTF stacks can also 
be seen in Fig. O which for each band state shows the 
electronic charge localized on the TCNQ and the TTF 
molecules, respectively. 

The metallic nature of the TTF-TCNQ band struc- 
ture arises from the energetic overlap of the quasi-lD 
TCNQ and TTF bands and the electronic charge trans- 
fer between them. Due to interstack interaction our cal- 
culation yields two Fermi vectors, kj^ = 0.27 and 0.33 
A~^, slightly larger than expected from the nesting vec- 
tor 2kF = 0.485 A-i derived from the CDW periodicity. 
The theoretical bandwidths along FZ are 0.7 eV (TCNQ) 
and 0.65 eV (TTF), in fair agreement with experimental 
estimates of ^0.5 eV.^ 

The conduction band dispersion perpendicular to b* 
is essentially negligible. Along the FY line (a*) it is 
practically zero. Slight dispersion of the TCNQ-derived 
band occurs along FB, which results from weak interac- 
tion along c between the molecular end groups in neigh- 
boring TCNQ stacks. However, the effect is too small 
to cause a band crossing along FB. The resulting Fermi 
surface topology is therefore truly ID. 



IV. ANGLE-RESOLVED PHOTOEMISSION 

ARPES measurements have been performed at our 
home lab, using He I radiation from a discharge lamp 
and an Omicron EA 125 HR electron energy analyzer, 
and with synchrotron radiation at BESSY (Berlin) us- 
ing an Omicron AR 65 spectrometer In both cases the 
energy and angular resolution amounted to 60 meV and 
±1°, respectively. All data were taken above the Peierls 
transition at a sample temperature of 60 K. TTF-TCNQ 
single crystals were grown by diffusion in pure acetonitrilc 
and had typical dimensions of 2 x 5 x 0.2 mm'^, with the 
long axis along b. Their quality was characterized by x- 
ray diffraction, electron spin resonance, and conductivity 
measurements. Clean surfaces parallel to the ab-plane 
were obtained by in situ cleavage of the crystals at a base 
pressure of < 10~^° mbar. The stoichiomctry of this sur- 
face, which contains both TCNQ and TTF chains,"!^ was 
verified by x-ray photoemissioni^ From the observation 
of momentum-dispersive ARPES structures we conclude 
on a crystalline long-range order of the surface, which 
has independently been confirmed by scanning tunneling 
microscopy (STM)i2Si2i 

Great care was taken to avoid photon-induced surface 
damage by minimizing the exposure to the incident 25 
eV radiation. The effect is demonstrated in Fig. 21 which 
contains spectra taken at the experimental Fermi vector. 
For a freshly prepared surface, i.e. immediately after 
cleavage of the crystal, an intense peak is observed close 
to the Fermi level. After two hours of exposure to vacuum 
ultraviolet (VUV) radiation its intensity has strongly de- 
creased and its peak position shifted by more than 0.1 



eV away from Ep. However, the original spectrum can 
be completely recovered by taking data on another pre- 
viously unexposed sample spot. This demonstrates that 
the observed surface degradation is not simply due to 
contamination or decomposition in the vacuum but in- 
deed caused by VUV radiation. Unfiltered higher order 
light or direct use of higher photon energies (^35 eV) re- 
duces the time scale of the VUV-induced surface damage 
even down to minutesiS All data presented in the remain- 
der of this paper were obtained before noticeable surface 
decomposition occured. 

Figure shows energy distribution curves obtained 
along the b* axis, i.e. the ID direction. The spectral 
features display pronounced dispersion, whereas spec- 
tra measured perpendicular to b* are dispersionless (not 
shown here, see Refs im . Our data are in excellent 
agreement with those of Zwick et al.^^ but display in 
parts more spectral detail. For example, in normal emis- 
sion {0 = 0°) we can clearly distinguish two peaks at 
0.19 and 0.54 eV below Ep, labeled a and b in Fig. 
In Ref. peak a appeared only as a shoulder and was 
not discussed. For off-normal emission peak b splits into 
two parts. The upper one (retaining the label b) moves 
upwards in energy and converges with a close to = 6°, 
where both features reach their closest approach to the 
Fermi level. We identify this position as Fermi vector 
which yields kp = 0.24 ±0.03 A~^, in good agreement 
with the value derived from the CDW vector. Note how- 
ever, that despite the high conductivity no metallic Fermi 
edge is observed in the fcp-spectra (see also Fig.QJ, as al- 
ready observed in Ref. [T^ The spectral intensity rather 
decreases almost linearly down to zero at exactly the 
Fermi level. 

Beyond kp a weak structure c moves back again from 
the Fermi level and displays a dispersion symmetric 
about 6 — 22° corresponding to the Z-point of the Bril- 
louin zone»2& Returning to the splitting of peak b away 
from 9 = 0° we note that its lower part (labeled d) dis- 
perses downwards in energy and eventually becomes ob- 
scured by peak c. For very high emission angles, corre- 
sponding to a k-vector in the next zone, one observes a 
symmetry-related weak shoulder d'. 



V. TEMPERATURE DEPENDENCE OF THE 
PHOTOEMISSION SPECTRA 

Upon cooling through the Peierls temperature the 
Fermi vector spectrum has been shown to display the ex- 
pected opening of a CDW gap.^^ However, even more 
remarkable is the temperature dependence of the kp- 
spectrum above the transition^Sii^ This is shown in 
Fig. inia) where considerable spectral changes are ob- 
served between 60 and 260 K. The interpretation of these 
changes hinges on a careful intensity normalization of the 
spectra. This has been achieved by normalizing them on 
the residual background intensity above the Fermi en- 
ergy, which is a good measure of the exciting photon 
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flux.^^ Incidentally, this procedure leads to an almost 
complete alignment of the temperature spectra at high 
binding energies (^1.3 eV). At lower energies the spec- 
tral changes from low to high temperatures can then be 
described as a pronounced intensity loss of the peak near 
the Fermi level (and its slight shift away from Ep), while 
at the same time the intensity increases between -0.4 and 
-1.3 eV (see difference spectra in Fig. Elb)). With the 
described normalization the integrated spectral weight 
remains however conserved within experimental uncer- 
tainty. We also note that these temperature effects are 
fully reversible. 

These observations and, in particular, the conservation 
of the total intensity suggest that At k = kp spectral 
weight is transferred from low to high binding energies 
with increasing temperature. However, based on the tem- 
perature dependence of the fci^^-spectrum alone we can- 
not rule out the possibility that the effect is caused by 
a redistribution of spectral weight in momentum space 
rather than in energy, caused, e.g., by phonon-induced 
disorder which may be large in organic compounds. In 
order to check this we have determined the /c-integrated 
density of states (DOS) by summing up ARPES spec- 
tra covering the entire ID Brillouin zone from F to Z. 
Fig. EIc) shows the resufr for T = 60 and 300 K. A 
smearing of the spectral weight distribution in momen- 
tum space due to thermally excited phonons should have 
no effect on the fc-integrated energy spectrum, except 
possibly for the phonon-induced lifetime broadening of 
the spectral peaks. The latter seems to be the case 
for the temperature change around —1.6 eV (the corre- 
sponding ARPES spectra show that the broadening oc- 
curs only near the zone edge, indicating a particularly 
strong electron-phonon coupling there). At low binding 
energies we recover the temperature dependence of the 
hp spectrum (cf. the difference spectra in Fig. Elb) and 
|Sfd)), which can clearly not be explained by line broad- 
ening. We hence conclude that the temperature behavior 
at the Fermi vector is indeed caused by a spectral weight 
transfer in energy. 



VI. COMPARISON OF PHOTOEMISSION AND 
BAND THEORY 

The identification of the ARPES dispersions indicated 
by the thin lines in Fig. |5l is further substantiated by 
a different representation of the data. Figure [7| shows 
the negative second energy derivative of the photocur- 
rent —(Pl/dE^, clipped at zero value, as grayscale plot in 
the {E, fc)-plane. This enhances the visibility of the spec- 
tral structures and visualizes their dispersion in a com- 
pletely unbiased way. Also shown are the DFT conduc- 
tion bands. The comparison of experiment and theory re- 
veals qualitative similarities but also significant discrep- 
ancies. Starting our discussion with experimental struc- 
ture c we find its dispersive behavior in agreement with 
that of the theoretical TTF-derived bands, except that 



the (occupied) bandwidth exceeds the theoretical one by 
a factor of ^ 1.7. Similarly, structures a and b can be 
attributed to the theoretical TCNQ doublet bands if one 
accounts for largely enhanced (a: ^ 2.0 and b: ^ 2.4) 
band widths. Finally, we point out that experimental 
feature d finds no counterpart in the band calculation. 

Our experimental conduction band widths are not only 
at variance with band theory but also clearly exceed the 
estimates derived from bulk-sensitive measurements]^ As 
the ARPES probing depth is comparable to the thickness 
of a single molecular layer (c/2 = 9.23 A, cf. Fig. IS)), the 
observed discrepancies suggest that the electronic struc- 
ture of the topmost layer(s) differs from that of the vol- 
ume. One possible origin could be a structural surface re- 
laxation involving the tilt angles of the planar (and rela- 
tively rigid) TTF and TCNQ molecules, respectively, rel- 
ative to the b-axis. We note that these angles correspond 
to a total energy minimum configuration resulting from 
a competition between maximum covalent bonding along 
the stack direction and minimum Coulomb energy in the 
Madelung potential of the surrounding molecular ionsi^ 
It seems conceivable that at the surface this balance is 
offset due the altered Madelung potential, leading to dif- 
ferent equilibrium tilt angles of the topmost molecules, 
as sketched in Fig. |H1 If the surface tilt happens to be 
larger than that in the volume, it will result in a reduced 
separation between the molecular planes within a given 
stack. This in turn leads to an increase of the intermolec- 
ular hopping integral and hence the bandwidth. 

We have tested this idea by performing a band calcu- 
lation for a hypothetical volume structure with increased 
tilt angles (approximately doubled relative to their bulk 
values), which indeed leads to a strongly enhanced band 
width at least for the TTF-derived bands in good agree- 
ment with their ARPES dispersion. However, a realistic 
calculation of the surface relaxation by total energy op- 
timization for a semi-infinite crystal is currently beyond 
our technical limits owing to the large size of the TTF- 
TCNQ unit cell. Unfortunately, a reliable experimen- 
tal determination of the molecular surface tilt seems also 
out of reach, as the usual methods for surface structure 
determination do not work here. Low energy electron 
diffraction (FEED) of TTF-TCNQ is strongly hampered 
by electron-induced surface damage even faster than that 
caused by the VUV photonsiS STM as another impor- 
tant structural surface probe is only capable of deter- 
mining the surface periodicitj^S*^ but cannot give any 
reliable information on molecular off-plane orientation. 
Therefore, the suggested surface relaxation has to re- 
main a mere speculation at this point. However, what- 
ever its microscopic origin, the observed enhancement of 
the ARPES bandwidth with respect to the volume is an 
experimental fact and we thus have to accept it as an 
established property of the probed surface layer which 
distinguishes it from the bulk. 

We are finally left with feature d, which even under 
the assumption of a surface band width renormalization 
cannot be identified with any of the theoretical volume 
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bands. It might appear tempting to attribute it to an 
intrinsic surface state. However, such an interpretation 
is in conflict with the observed Fermi vector of the other 
bands: Since d stays wefl below the Fermi level and would 
thus be occupied throughout the entire Brillouin zone 
(in a one-electron band picture), it should severely affect 
the delicate charge balance between the TTF and TCNQ 
bands and shift the surface Fermi vector notably from its 
bulk value, which is not the case. Another explanation 
of d as backfolded image of the TTF band induced by 
long-ranged CDW fluctuations^i is ruled out due to the 
lack of other evidence for backfolding in the data. As 
we will discuss in the following section, spectral feature 
d finds a natural explanation as a many-body effect. 



VII. COMPARISON TO THE ID HUBBARD 
MODEL 

There is substantial experimental evidence that 
Coulomb interaction plays an essential role in the elec- 
tronic structure of TTF-TCNQL2iSi2^ and that a purely 
band theoretical description may be inadequate. On the 
theoretical side, the dramatic effects of electron-electron 
interaction on the low-energy properties of ID metals 
have been studied in much detail using the Tomonaga- 
Luttinger (TL) model. ^ It is based on a ID conduction 
band with infinite linear dispersion and treats interac- 
tion by including scattering processes about the Fermi 
pointsii*^ The TL model focusses on the low-energy 
physics and describes in detail the breakdown of the 
quasiparticle picture for the low-lying excitations and the 
emergence of TLL behavior resulting from the dynami- 
cal decoupling of spin and charge. For example, the TL 
model predicts a low-energy onset of the single-particle 
spectrum which is no longer given by a metallic Fermi 
edge but rather by a power law behavior oc oj" , with the 
exponent a determined by the coupling parameters of the 
model. The low-energy physics of the TL model defines in 
fact a universality class which includes also more compli- 
cated ID models of interacting electrons i^^i^ However, 
by its very construction the TL model contains no in- 
trinsic energy scale, and therefore the energy range of its 
applicability to real ID metals is principally unknown, 
making it less useful for the study of finite-energy spec- 
tral properties. 

As seen in the previous section, the ARPES data of 
TTF-TCNQ indeed show unusual behavior on an en- 
ergy scale of the entire band width, not just near the 
Fermi level. The spectral properties over this much 
wider energy range have so far only been addressed 
by the ID single-band Hubbard model. Compared to 
band theory it appears as a much better starting point 
for the description of TTF-TCNQ and other organic 
charge transfer salts. In fact, various properties of these 
ID conductors have already successfully been analysed 
within a Hubbard model framework, such as the magnetic 
susceptibilitySMZ or the nuclear spin relaxation ratei2i2£ 



The underlying idea is that the local interaction energy 
U for two electrons residing on the same molecule domi- 
nates over long range Coulomb contributions. The de- 
localization of the charge carriers is described by the 
hopping integral t (the bare bandwidth amounts to At 
in one dimension). The Hubbard model also defines 
an intrinsic energy scale for spin excitations, which for 
large values of U/t is given by the exchange constant 



2f^ 
U 



J = 

parameter [n 



i(27rn 



2tt 



with n being the band filling 

0.59 for TTF-TCNQ). 

At low excitation energies the physics of the ID Hub- 
bard model with finite doping follows the TLL phe- 
nomenology. The C//i-dependence of the non-classic TLL 
exponents, which control the asymptotics of the low- 
energy correlation functions, can be extracted from its 
Bethe-ansatz solutioniii However, in contrast to the TL 
model the Hubbard model also allows the study of finite- 
energy properties. Recently, an exact analysis on the 
basis of the Bethe-ansatz has showr>iS*ii that all energy 
eigenstates of the ID Hubbard model can be described 
in terms of occupancy configurations of various collective 
spin-only and charge-only modes, namely spinons (zero- 
charge spin excitations), holons (spinless charge excita- 
tions), and a third type of charged quantum objects^S 
We refer the reader to Refs. "T^andl?' for details. The im- 
portant point is that this description is valid for all energy 
scales of the model and follows from the non-perturbative 
organization of the electronic degrees of freedom. 

Here we are interested in the electron removal spec- 
trum of the ID Hubbard model. Qualitative properties 
of the spectrum have already been derived from early 
calculations within the strong coupling limit {U/t — > 
QQ-jj43j44 ]y[Qj.g recently, it has become possible to deter- 
mine the spectral behavior also for intermediate interac- 
tion strengths {U ^ 4i), either by numerical methods^ 
or by Bethe-ansatziiSiii Here we will restrict ourselves 
to the energy- vs.-momentum dependence of the spec- 
tral features, which can be calculated exactly with the 
latter method. The calculation of matrix elements be- 
tween ground and excited states and hence of the spec- 
tral weight distribution is more complicated with this 
method and will be presented elsewhere. A schematic 
picture of the spectral dispersions is given in Fig. I^a). 
As the hole generated by the removal of one electron de- 
composes (or "fractionalizes" ) into decoupled spin and 
charge excitations, there is a manifold of ways to dis- 
tribute excitation energy and momentum among these 
collective modes giving rise to an excitation continuum, 
indicated by the shaded area in Fig. Ufa). 

However, due to the phase space available for electronic 
hole fractionalization this continuum is not structureless. 
It is dominated by lines of singularities (solid curves in 
Fig. inia)) which roughly speaking correspond to situa- 
tions, in which either the charge mode propagates with 
the entire excitation energy leaving zero energy for the 
spin channel, or vice versa. We denote these dispersion 
curves hence as "charge" and "spin" branches, respec- 
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tively. At the Fermi vector both branches are degener- 
ate, but due to their different group velocities they spht 
away from the Fermi level. This low-energy behavior 
has already been found for the TL model. ^ The ID Hub- 
bard model now allows us to explore also the finite-energy 
dispersion of these features. The spin branch for exam- 
ple reaches its maximum binding energy for momentum 
A: = at about ^J, reflecting the dispersion of a bare 
spinon. The dispersion of the charge part is a little bit 
complicated. Starting from k = —kp it follows a nearly 
parabolic-like dispersion reaching at fc = ~\-kp a maxi- 
mum binding energy which scales with the hopping in- 
tegral t. From there it disperses upwards again under a 
pronounced loss of spectral weighl4ii4i^ until it even- 
tually crosses the Fermi level at -\-Zkp (for symmetry 
reasons there is a corresponding charge branch running 
from +kp to —ikp). The shape of its dispersion reflects 
that of a bare holon, with the distance of Akp (rather 
than 2kp) between its Fermi level crossing points owing 
to the fact that the holon is a spinless quantum object. 
The shift of the symmetry point away from fc = to ±kp 
can be understood from a detailed microscopic analysis 
of the electronic hole spectrumiiSiii The peculiar high- 
energy behavior of the charge branch was first noted by 
Penc et al. in the case U /t oo.^'^ 

Comparing this picture to the observed experimental 
dispersions of the TCNQ-related peaks in the ARPES 
data (cf. Fig.!?)) we find remarkable similarities. In fact, 
it is even possible to obtain a quantitative Hubbard model 
description of the experimental dispersions. For this pur- 
pose we have utilized the Bethe-ansatz method intro- 
duced in Ref. [461 Further details about the line shape 
predicted by the Hubbard model within such a method 
for the TCNQ dispersions are presented elsewheren^ 
The method leads to U /t dependent branch lines which 
are given by the expressions of Ref. il6( [Eqs. (C15), 
(C16), (C19), and (C21)] for the bare holon and spinon 
dispersions,™ which in turn reproduce those of the charge 
and spin branches in the electron removal spectrum. For 
these calculations the model parameters U and t were 
chosen in such a way as to yield optimum agreement with 
the ARPES dispersions. The comparison of the model 
calculation to the experimental TCNQ dispersions (from 
ARPES spectra measured on a finer fc-grid than those 
presented above) in Fig. Elb) yields an almost perfect 
match. This allows us to identify the experimental struc- 
tures a and b (cf. Fig. d) as spin branch and the upper 
part of the charge branch, respectivelyi^ Moreover, the 
as yet unidentified structure d now finds its natural expla- 
nation as the high-energy part of the theoretical charge 
branch, at least for not too large fc-vectors. Experimen- 
tally, its reversed dispersion beyond kp and its eventual 
3fci?-crossing is not observed, most likely due to the the- 
oretically predicted loss of weight at larger fc and the 
overlapping TTF band. The model parameters used to 
fit the theoretical dispersions to the experimental ones 
are t = 0.4 eV and U — 1.96 eV, corresponding to a 
rather moderate coupling strength oiU /t = 4.9. 



We finally turn to the TTF-related ARPES feature. 
Concerning its dispersion we observe no extraordinary 
behavior other than the enhanced bandwidth relative to 
band theory. Complimentary to the TCNQ band, the 
TTF-derived conduction band (or rather band doublet) 
is more than half-filled (n = 2 — 0.59>1). For this case 
the Hubbard model predicts a charge branch line whose 
dispersion shows some similarities to that of the exper- 
imental feature c in Fig. [3 However, there should also 
be an additional spin branch line for which we observe 
no clear evidence. It is possible that the model parame- 
ters suitable for the TTF band differ from those for the 
TCNQ band and are such that spin-charge separation is 
less pronounced. Moreover, x-ray diffraction studies have 
reported the existence of Akp CDW fluctuations on the 
TTF chains up to 220 K.®'"'^ It is not clear how such 
fluctuations may affect the ARPES spectra; they could 
for example account for the relatively large line-width 
of the experimental TTF peak (cf. structure c in Fig.[SJ), 
thereby obscuring a possible small spinon-holon splitting. 
At this point the interpretation of the TTF part of the 
ARPES data has to remain an open question. 



VIII. DISCUSSION 

The ID Hubbard model thus provides a quantitative 
description of the experimental TCNQ-related disper- 
sions and an explanation for the failure of band the- 
ory. In fact, earlier studies of the electronic and mag- 
netic low-energy properties of TTF-TCNQiSiSSiS have 
already used this model successfully for the interpretation 
of their data. They estimated that the local interaction 
energy U and the bandwidth At are comparable and of 
the order of 1 eV. This is consistent with the parameters 
of our model fit. Concerning the resulting bandwidth 
At = 1.6 eV we observe an approximate doubling with 
respect to the result of our DFT calculation (0.7 eV), 
just as in the bare band-theoretical interpretation of the 
ARPES data. It is again attributed to a possible molecu- 
lar surface relaxation as already discussed in Section IVIl 
The hopping integral of our Hubbard model fit thus re- 
fiects a surface property. In order to compare to bulk 
properties we should rather use the value inferred from 
the DFT bandwidth, t = 0.175 eV. As the intramolec- 
ular Coulomb energy [/ is a local quantity, we do not 
expect large differences between bulk and surface. With 
U = 1.96 eV we thus obtain for the coupling strength 
in the volume a value of U/t = 11.2. From U and t we 
can also calculate the magnetic exchange energy J, which 
for the volume yields 21 meV (110 meV for the surface). 
This is in good agreement with experimental estimates 
of J inferred from magnetic susceptibility measurements 
which range between 17 and 30 meVi22iS. 

Further evidence for the importance of correlation ef- 
fects is provided by the unusual temperature dependence 
of the photoemission spectra. Commonly, temperature 
effects are caused by electron-phonon interaction with 
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spectral changes occuring on an energy scale ksT, due 
to an altered population of phonons with a comparable 
energy (at least within harmonic approximation and with 
linear coupling). However, for TTF-TCNQ we observe 
upon warming-up a shift of spectral weight from low to 
high binding energies by ~ 1 eV, i.e. an energy of the 
order of the bare bandwidth and hence much larger than 
the thermal energy scale. This seems to rule out con- 
ventional electron-phonon coupling as the origin of the 
temperature dependence, though we cannot completely 
exclude the possibility of non-linear coupling effectsi^ 

A much more natural explanation of the observed tem- 
perature effects can be inferred from calculations for the 
quarter-filled ID tJ model in the strong-coupling limit 
{J/t 0, corresponding to the U/t — ^ oo case of the 
Hubbard model) Here it was found that, compared 
to the zero temperature spectrum, considerable spectral 
weight is redistributed from the " spin" peak at the Fermi 
level to the bottom of the "charge" band at —2t already 
at temperatures < fc^T << t, exactly as observed in 
our data. We are not aware of similar calculations for 
moderate interaction strengths, but we expect this result 
to hold qualitatively also for finite J/ 1 or U/ 1, respec- 
tively. 

In conclusion, both the dispersive behavior of the 
TCNQ-derived ARPES structures as well as the tem- 
perature dependence of the spectra are found to be well 
accounted for, in parts even quantitatively, by the finite- 
energy spectral properties of the ID Hubbard model. The 
observed discrepancies to band theory thus appear as 
a consequence of spin-charge separation, which occurs 
in that model on all energy scales. In this interpreta- 
tion our ARPES resuhs on TTF-TCNQ represent the 
first spectroscopic observation of spin-charge separation 
in a quasi- ID metal on an energy scale of the conduction 
band width. It is interesting to note that there exists 
independent experimental support for the occurence of 
spin-charge decoupling in TTF-TCNQ from the contrast- 
ing temperature dependence of conductivity and spin 
susceptibilityi^ 

We close this section with a discussion of the spectral 
onset near Ep, for which the ID Hubbard model predicts 
a low energy behavior cx |£' — Epl" with the exponent 
ranging between a = and a = 1/8 for U/t ^ and 
U/t — > CX), respectively^ This is in clear contrast to our 
experimental observation of an almost linear energy de- 
pendence, for which there are various possible explana- 
tions. First of all, the physics of the ID Hubbard model 
is expected to be applicable only for excitation energies 
larger than the transverse transfer integrals associated 
with interchain hopping (the DFT band dispersions of 
Fig. El give an estimate of the relevant energy scale). In 
addition, long-range interactions beyond a simple Hub- 
bard model (e.g. the effect of nearest neighbor interac- 
tion, which may be non-negligible for TTF-TCNQ&2.) are 
capable to increase the onset exponent up to a T— Un- 
fortunately, the spectral properties of extended Hubbard 
models at higher binding energies are not well known. Fi- 



nally, it has recently been argued that low-energy power 
law exponents of the order of unity can also be caused 
by impurities and/or defects on the surface of an organic 
conductor, which localize the ID electrons to strands of 
finite length, leading to the concept of a "bounded Lut- 
tinger liquid" 

However, it seems likely that the failure of the simple 
Hubbard model at low energies is not just a purely elec- 
tronic effect. Rather, one should also expect pronounced 
contributions by electron-phonon coupling, which after 
all is strong enough to drive a Peierls transition at low 
temperatures. On the other hand, a simple interpretation 
of our linear spectral onset in terms of a Peierls pseudo- 
gap due to CDW fluctuations above Tp = 54 K must be 
ruled out, as the size of the underlying low-temperature 
Peierls half-gap is only 20 meV,'* much smaller than 
the energy range of the onset (~ 100 meV). Further- 
more, 2kp CDW fluctuations disappear already at 150 
K, while Akp fluctuations - though still observable at 300 
K - strongly weaken with increasing temperature>^ In 
contrast, the ARPES spectral weight near Ep is found 
to become reduced (while still being linearly energy- 
dependent) from low to high temperatures (cf. Fig. 
The large energy range of the spectral onset indicates the 
importance of coupling to other phonons than those in- 
volved in the Peierls transition and is consistent with the 
phonon spectrum of TTF-TCNQ, which indeed reaches 
up to - 200 meV.^^ Even so, any detailed understanding 
of the spectral properties of TTF-TCNQ at low energies 
will require the consideration of electronic correlations 
and electron-phonon coupling effects on an equal foot- 
ing, which remains to be a challenge to modern solid state 
theory. Whatever the details of any such description, our 
above results indicate that its high-energy physics must 
be close to that of the ID Hubbard model. 



IX. CONCLUSION 

The electronic structure of TTF-TCNQ above the 
Peierls transition as probed by ARPES deviates signif- 
icantly from DFT band calculations. The experimen- 
tal observation of an approximate doubling of the over- 
all conduction band width relative to band theory is at- 
tributed to a structural relaxation of the topmost molec- 
ular layers. When accounted for an enhanced electron 
hopping integral at the surface, the spectra of the TCNQ- 
derived bands can be brought into consistent and even 
quantitative agreement with the theoretical flnite-energy 
single-particle spectrum of the ID Hubbard model. This 
picture is further supported by a temperature-dependent 
redistribution of spectral weight over energies much 
larger than the thermal energy. Within this interpreta- 
tion our experimental results provide spectroscopic evi- 
dence for spin-charge separation on an energy scale of the 
conduction band width. In contrast, the spectral behav- 
ior at low binding energies is found to deviate from that 
of the simple ID Hubbard model, possibly due to higher 



8 



dimensional effects combined with the additional impor- 
tance of electron-phonon coupling and possibly also long- 
range electron-electron interaction. TTF-TCNQ thus 
represents an interesting model system to study elec- 
tronic correlation effects in a ID metal. 
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FIG. 1: Crystal structure of TTF-TCNQ. ■dp and i?q indicate the tilt angles of the planar TTF and TCNQ molecules, 
respectively, relative to the ac-plane. Also shown is the monoclinic Brillouin zone with its high symmetry points. 
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FIG. 2: DFT band structure near the Fermi level along the three major high-symmetry lines of the Brillouin zone of TTF-TCNQ. 



12 



TCNQ bands 



LU 

O -1 h 



]88888888880OQoooo 








-3 - 



oooA 



^nooooooooi 



,0 



2a °° 



ooo 



z r 



FIG. 3: Theoretical band dispersions along FZ showing the molecular origin of the bands. The size of the symbols represents 
the charge of each state residing on the TCNQ (left panel) and TTF (right panel) molecules. 



13 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




energy relative to (eV) 



FIG. 4: Angle-resolved energy distribution curves at the Fermi vector showing the effect of photon-induced surface degradation 
{hu = 25 eV, T = 61 K). For a detailed discussion see text. 
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FIG. 5: Energy distribution curves measured along the FZ direction (hv = 25 eV, T = 61 K). The thin lines are guides to the 
eye and are meant to indicate the dispersion of the spectral features. 
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FIG. 6: Temperature dependence of the photoemission spectra {hv = 21.2 eV): (a) Momentum-resolved spectrum at = feF 
measured between 60 K (blue/dashed curve) and 260 K (red/solid curve), (b) Difference spectra relative to 60 K. (c) Momentum- 
integrated spectrum at 60 K (blue/dashed) and 300 K (red/solid); note the larger energy scale compared to (a) and (b). (d) 
Difference spectrum relative to 60 K. 
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FIG. 7: Gray-scale plot of the ARPES dispersions (see text for details) in comparison to the conduction band dispersions 
obtained by density functional theory. 
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FIG. 8: Schematic picture of a possible surface relaxation leading to enhanced molecular tilt angles i?q,s and i?f,s of the topmost 
TCNQ and TTF molecules, respectively. 
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FIG. 9: (a) Schematic electron removal spectrum of the doped ID Hubbard model with band filling 1/2 < n < 2/3. The shaded 
region denotes the continuum resulting from spin-charge separation. The solid curves indicate the dispersions of the "spin" 
and "charge" singularity branches (see text for details), (b) Theoretical spin and charge branch dispersions of the ID Hubbard 
model calculated for U = 1.96 eV, t = 0.4 eV, and n — 0.59 in comparison to the ARPES dispersions of the TCNQ-derived 
conduction band complex (measured with He I radiation, representation as in Fig. 0. 



